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TECHNICAL  MMORANDUM  1260 


EXACT  SOLUTIONS  OF  EQUATIONS  OF  GAS  DYNAMICS* 

By  I.  A.  Kiebel 

The  equations  of  the  tvo -dimensional  stationary  problem  of  gas 
dynamics  are  of  the  form 
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where  u  and  v  are  the  components  of  the  velocity  along  the 
coordinate  axes  x  and  y,  respectively,  p  is  the  pressure, 
p  the  density,  and   K=  c_/c^,  the  ratio  of  specific  heats. 
The  equation  of  continuity  permits  the  introduction  of  the  stream 
function  from  the  equations 
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The  last  of  equations  (1)  gives 
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♦"Primer  Tochnogo  Reshenia  Ploskoi  Vikhrevoi  Zad.achi  Gazovoi 
Dinamiki."     Prikladnaya  Matematika  i  Mekhanica.     Vol.  XI,    1947, 
pp.   193-198. 
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where  ^  is  a  certain  function  only  of  ^.     The  first  two  of 
eq,uation3  (1)  give  the  Bernoulli  law 


(4) 


where  i^  is  a  function  of  ^^     and  the  equation  for  the  vorticity 
Is 

dx  ■  Sy  ■  ■  ^\d^i'  "k-1  ^    dv|//  ^^^ 


For  the  solution  of  the  vortex  problem  in  which  at  least  one 
of  the  derivatives  diQ/d>i/  and  d^/d'4'  is  different  from  zero,  it 
is  convenient  to  pass  from  the  variables  x  and  y  to  the  vari- 
ables T*  =  X     £Lnd  ^|/.  Eq.uation3  (2)  and  (5)  then  assume  the  form 
(See,  for  example,  reference  1.) 
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The  problem  reduces  to  the  determination  of  the  five  functions 
u,  V,  p,  p,  and  y  of  x*  and  '^     from  eq.uation3  (3),  (4), 
and  (6). 


By  using  the  last  of  equations  (6),  the  function  X  is  intro- 
duced with  the  aid  of  the  equations 
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(7) 
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The  first  two  equations  of  equation  (6),  on  the  "basis  of  equa- 
tions (3)  and  (4),  then  assume  the  forms 
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where  the  asterisk  on  x  has  been  dropped. 

By  differentiating  equation  (8)  with  respect  to  ^  and  equa- 
tion (9)  with  respect  to  x,  j     can  be  eliminated  and  a  single 
equation  for  the  function  X  obtained. 


In  order  to  obtain  an  example  of  the  exact  solution  of  the 
system  of  equations,   i^  is  set  equal  to  a  constant: 


X  =  -  H  (\^)  X 


K-1 
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(10) 


where  H  is  a  certain  function  of  'I'  to  be  determined.  Equa- 
tions (8)  and  (9)  then  give 
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When  equation  (ll)  is  differentiated  with  respect  to  '^      and 
equation  (12)  is  differentiated  with  respect  to  x,  terms  with 
the  same  degree  of  x  are  collected  and  after  simple  trans - 
formations  there  is  obtained 

_1 

Wl     J  2  ^       (13) 

H'  (H'^+KHH")'  =  -  (k-1)  H"  (H'^+KHH")^ 
Successive  integration  of  the  second  of  equations  (13)  gives 


H'^  +  KH  H"  =  CiH'-'-"'^ 


J5_   /  (14) 

H  =  Co  (C,-H'^+'*' 


'2^-1 

From  the  last  equation,  the  relation  between  "^    and  H  can 
"be  found  with  the  aid  of  quadratures.   It  is  more  convenient,  how- 
ever, to  introduce  a  numbering  of  the  streamlines  directly  with 
the  aid  of  H'   and  not  '^. 

Thus,  equations  (14)  give  H  as  a  function  of  H'   and,  using 
equations  (13),  permit  finding   ^  in  terms  of  H'   in  the  form 

l-K  2n 

Substituting  in  equation  (ll)  and  integrating  yields,  for  the 
streamlines. 


K-1  _2  »«^  _1 

y  =  I  H-x*^""^  [21^  -  CiH'^-'^x  "^^^     ^   dx  +  F  (E-)     (16) 
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Comparison  with  equation  (12)  shows  that  F  =  constant  and 
without  loss  of  generality  can  he  taken  as  F  s  0. 

For  the  determination  of  u,  v,  p,  and  p  the  following 
relations  are  written : 


u^  =  21^  -  C]^E'^"'^x 
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For  simplification  it  is  convenient  to  introduce  the  nondimen- 
sional  variables  x  and  y  and  in  place  of  H'   introduce  t\ 
with  the  aid  of  the  following  equations; 


X  =  Lx 
y  =  Ly 


(20) 
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where 
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Equations  (16)  ajid  (17)  then  assume  the  following  simple  form: 
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gives 


If      K  =  1.4,     the  integral  in  equations   (22)    is  evaluated  and 
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(24) 


and  the  streamlines  can  therefore  he  easily  constructed. 

Figure  1  shows  the  streaml-nes  for  t^  ±  1,  2,  3,  .  .  .  ,  19. 

This  motion  possesses  hoth  supersonic  and  subsonic  velocities, 
for  the  line  of  transition  (shown  dotted  in  fig.  1)  is  obtained  if 
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It  is  possible,  without  difficulty,  to  construct  the  character- 
istics in  the  x,y-  and  u,v-planes.   Instead  of  the  equation  of 
the  epicycloid  in  the  u,v-plane,   in  the  case  considered  a  more 
complicated  equation  arises. 

In  the  vortlclty  problem,  along  the  characteristic  there 
occurs 

ctg  a       sin  a  cos  a       ,  ,   , 

dp  T  — - —  dw  =  ±  d  log  ^  (26) 

w  K-1 

vhere  a  is  the  Mach  angle,  w  the  magnitude  of  the  velocity,  and 
P  the  angle  of  Inclination  of  the  velocity  to  the  x-axis.  Equa- 
tions (23)  yield 

_   21^-u^    21.-w^  cos^  p  ■ 

n     =   °     =  —2 (27) 

1^       v^  sln^  p 

But  ^   depends  only  on  H',  that  is,  on  t^  .  According  to 
equation  (15),  the  following  is  obtained: 


^  =  1^  (=^  #""=1  T  <^^-'^ 


,       .ill      p  K-1         2K 
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Substituting  this  value  of  ■d   in  the  right  side  of  equa- 
tion (26)  yields  the  differential  equation  for  the  characteristics 
In  the  PjW-plane. 

Finally,  the  question  arises  whether  it  is  possible  to  obtain 
such  vortex  motion  by  transition  through  a  surface  of  strong  dis- 
continuity of  some  other  kind  but  with  Irrotatlonal  motion.  This 
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problem  may  "be   answered  in  the  affirmative.  For  on  a  surface  of 
discontinuity  there  must^  among  others,  be  satisfied  the  relations 
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where 
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"+j 


and  0^  are  the  pressure,  density,  and  velocity 


of  propagation  of  the  surface  of  discontinuity  on  one  side  of  this 
surface  and  p_  emd  p_  are  the  pressure  and  density  on  the  other 
side.  The  magnitudes  p^  and  p^.  may  be  taken  from  the  vortex 
motion.  The  magnitude  6^     is  found  in  terms  of  the  elements  des- 
cribing the  motion  and  in  terms  of  the  inclination  of  the  surface 
of  discontinuity.  Finally,  p_  and  p_  can  be  connected  by  the 


relation  p_    =  ■diP_, 


where 


is  a  constant  (up  to  a  certain 


degree  of  arbitrariness)  magnitude.   (At  the  left  of  the  surface  of 
discontinuity  the  motion  is  irrotatlonal. ) 
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where  5  is  the  angle  between  the  normal  to  the  surface  of  dis- 
continuity and  the  x-axis  (the  normal  is  directed  toward  the 
"positive"  region);  M  may  be  expressed  in  terms  of  known  magni- 
tudes and  6. 

By  expressing  tan  6  in  terms  of  the  derivative  of  x 
with  respect  to  t|  along  the  surface  of  discontinuity  and  using 
equations  (22),  (18),  and  (19),  the  differential  equation  for 
determining  the  surface  of  discontinuity  (K  =  1.4)   is  obtained 
after  simple  transformations: 
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where  M  is  expressed  in  terms  of  i^   from  the  transcendental 
equation  (31)  in  which  the  constant  G  is,  to  a  great  extent,  arbi- 
trary. The  velocity  in  the  "negative"  region  will  be  determined 
first  on  the  surface  of  discontinuity  and  then  extended  on  the 
negative  region  by  the  usual  graphical  method  of  Busemann,  A  model 
of  the  motion  about  a  contour  having  an  angle  is  obtained. 

The  surface  of  discontinuity  extends  out  from  the  angle  and 
on  passing  through  the  angle  the  motion  reconverts  to  the  rotational 
motion  herein  considered. 

A  solution  analogous  to  that  developed  can  also  be  obtained 
for  the  problem  with  axial  symmetry.  Taking  for  the  independent 
variable  the  distance  r*  =  r  from  the  axis  of  symmetry  z  and 
the  stream  function  \|/  yields  the  relations 


br* 
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r*pvi,  J 


where  v^,  and  v,  are  the  velocity  components.  As  before, 
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Bernoulli's  lav  will  have  the  form 
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1/2   2\    Ki 


=  i. 


The  equations  of  Euler  give 

dr*  ~  ~SY 


The  function  X(r*,\l')  can  therefore  be  introduced  from  the 
conditions 
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Equations  (34)  and  (35)  now  permit  writing  (the  asterisks  on  r 
are  dropped) 
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where 
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By  eliminating  z  from  equations  (36),  a  single  equation  for 
the  determination  of  the  function  X  is  ohtained.   Particular 
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solutions,   analogous  to  the  solutions   in  the  first  case,    can  he 
constructed  hy  seeking      X   in  the  form 

X=   -  r      '^"'^  E(^) 


Translated  by  S.  Reiss 
National  Advisory  Committee 
for  Aeronautics . 
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